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Abstract 

A conservation law is one of the most fundamental properties in nature, but a certain class of 
conservation “laws” could be spoiled by intrinsic quantum mechanical effects, so-called quantum 
anomalies. Profound properties of the anomalies have deepened our understanding in quantum 
many body systems. Here, we investigate quantum anomaly effects in quantum phase transitions 
between competing orders and striking consequences of their presence. We explicitly calculate 
topological nature of anomalies of non-linear sigma models (NLSMs) with the Wess-Zumino-Witten 
(WZW) terms. The non-perturbative nature is directly related with the’t Hooft anomaly match¬ 
ing condition : anomalies are conserved in renormalization group flow. By applying the matching 
condition, we show massless excitations are enforced by the anomalies in a whole phase diagram in 
sharp contrast to the case of the Landau-Ginzburg-Wilson theory which only has massive excita¬ 
tions in symmetric phases. Furthermore, we find non-perturbative criteria to characterize quantum 
phase transitions between competing orders. For example, in 4 D, we show the two competing or¬ 
der parameter theories, CP( 1) and the NLSM with WZW, describe different universality class. 
Physical realizations and experimental implication of the anomalies are also discussed. 
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Introduction 

Quantum anomaly is one of the most fascinating phenomena in quantum many body 
systems. Quantum fluctuations spoil classical symmetry, thus corresponding conservation 
laws and Ward identities no longer hold and must be modified. Topological source terms to 
the “conservation” laws are induced by quantum fluctuations from anomalies, thus topology 
and symmetry are intrinsically tied. Consequences of the anomalies were first confirmed in 
the pion-decay (7r° —y 77), and since then deeper understanding has been achieved. 1 2 

Topological protection is one of the most fascinating properties of quantum anomalies, 
and remarkably this protection is independent of interaction strength, ’t Hooft first realized 
and applied these properties to confinement physics of quantum chromodynamics (QCD), 
so-called ’t Hooft mathcing, and constrained candidates of low energy degrees of freedom 
including the Goldstone bosons from the chiral symmetry breaking! 3 Such non-perturbative 
nature has been extensively applied to high energy physics, for example, the standard model, 
the Skyrme model of hadrons, and black hole physics.® 3 In condensed matter systems, it 
is also applied to several topological phases.® 19 For example, the presence of edge states in 
quantum Hall states and violation of the chiral current conservation in Weyl semimetals are 
examples of the realizations of the chiral 17(1) anomaly. Massless excitation in either edge 
or bulk is protected by anomalies’ topological nature. 

In this paper, we consider another realization of quantum anomalies in condensed matter 
systems, non-abelian anomalies in quantum phase transitions between competing orders 
and investigate consequences of their presence. We first show that a class of competing 
order theories has anomalies. Their presence becomes criteria to characterize competing 
order theories. With the criteria, we find that the competing order theory with the CP(1) 
deconhned transition in three spatial dimensions cannot describe the same universality 
class of the NLSM with WZW in sharp contrast to the case of two spatial dimensions where 
the two models are proposed to describe the same universality class. Furthermore, by 
using the ’t Hooft anomaly matching condition, we also find competing order physics with 
anomalies must contain massless excitation in sharp contrast to the case of the conventional 
Landau-Ginzburg-Wilson theory. We provide possible candidate theories of the quantum 
phase transitions with anomalies. 

Theories of Competing Orders 
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Various order parameters appear in strongly correlated systems, and intriguing interplay 
between order parameters has been reported. 20 2 ' In this section, we introduce two types of 
competing order theories, which describes fundamentally different competing order mecha¬ 
nisms, and set up our notation for later discussion. It is worthwhile to mention that we only 
focus on competing order theories in non-metallic systems in this paper. 

The phenomenological yd theory, so-called the Landau-Ginzburg-Wilson (LGW) theory, 
is one of the simplest ways to describe competing orders. Each order parameter (<pn, or tp M ) 
is a representation of corresponding symmetry group (say, vector representations of O(N) 
or O(M)). The minimal Landau functional is 

•Fl = r N if 2 N + u N <p% + r M cp 2 M + u m <Pm + s <PnVmi 

omitting fluctuation terms and higher order terms. Four different phases are basically de¬ 
scribed by the signs of the tuning parameters (r N , r M ) around the multi-critical point (0, 0) 
as shown in FIG. [lj Note that low energy excitation of the LGW theory’s symmetric 
ground state (S) is massive which can be easily shown by restoring fluctuation terms. One 
useful way to understand the phase diagram is to promote the symmetry O(N) x O(M) to 
0(N + M) symmetry by introducing a ‘super-spin’ (c PniVm ) and focus on the multi-critical 
point. Then, the four phases are accessed by introducing anisotropy operators. It is well- 
known based on symmetry that the theory is equivalently described by the non-linear 
sigma model (NLSM)J 2 ® 

„ , N+M 

S 0 = I d D x—(d(j) i ) 2 , 

g characterizes strength of fluctuations, and the four phases are again accessed by anisotropy 
operators near the critical coupling constant (g = g c ). 

Another type of competing orders theories can be obtained by incorporating topological 
terms to NLSMs. It is because NLSMs can easily include topological natures of the compact 
ground state manifold (S N+M ~ 1 ). There are two types of topological terms in NLSMs : 
NLSMs with the Wess-Zumino-Witten (WZW) term and NLSMs with the 0 term. The 
former is realized in D — N + M — 2, which is our main focus in this paper, and the latter 
is realized in D = N + M — 1, which will be discussed in future work. 
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FIG. 1: (a) Schematic phase diagram of two order parameters in LGW. The horizontal and 

vertical axes are tuning parameters of the order parameters. Four phases are determined by the 
signs of (rjv, tm) '■ S phase ({<Pn) = {'Pm) = 0), N phase {{<Pn) / 0, (pm) = 0), M phase 
((<Pn) = 0, (pm) A 0), and N + M phase ((pn) A 0, (<pm) A 0). The multi-critical point (0,0) 
is well described by a theory with an enlarged symmetry group (0(N + M )). For illustration, 
we choose N = 3 and M = 3 (for example, magnetism and valence bond solid) in three spatial 
dimensions, (b) Two generic phase diagrams with one parameter. The vertical axis is for magnitude 
of order parameters. The upper (lower) one is realized with the condition r n + > 0 (< 0), and 

the fine-tuned condition (rjy + = 0) gives a second order phase transition between the two 

symmetry broken phases. 

The 0(D + 2) NLSMs with the WZW terms in D space-time dimensions are 

S= f ( ^£+kS wzw , = (1) 

J.Xd 9 1=1 

with the spacetime manifold X D . The ground state manifold is S D+1 . The WZW term is 

r p. . p^i-'^D+i 

Swzw = ic D+ 1 / - ^d ^ i2 • • • , 

and the presence of the anti-symmetric tensor is related to the winding number of the ground 
state manifold, 7r£) + i(S' z?+1 ) = Z. We mainly focus on D — 2,4 though its generalization 
is straightforward. The boundary of the manifold Xr> + i is dXr> + i — Xj). The numerical 
constants are c D+1 = Area ^ g+l) and Area(S D+1 ) = ^+2)12) • ( c 3 = l and c 5 = ^). With 
a non-zero k, this class of competing order theories have been suggested to describe exotic 
competing order physics, for example, competition between Neel and valence bond solid 
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orders associated with deconfined quantum criticality] 29 30 

The engineering dimension of the coupling constant is [g 2 ] = 2 —D. In D > 2, the coupling 
constant is irrelevant indicating the model is well-suited to describe a weak coupling fixed 
point, a symmetry broken phase. Therefore, to access symmetric ground states, one need to 
consider strong coupling limits where a perturbative calculation is not reliable. Below, we 
show that non-perturbative nature of quantum anomalies allows us to investigate symmetric 
phases even in strong coupling limits. 

In D — 2, the coupling constant is marginal at the tree level, and Witten showed the 
WZW-NLSM is mapped to a massless fermion model by non-abelian bosonization. 31 This 
clearly shows that the WZW-NLSM in 2 D describes a different universality class from the 
LGW ip A theory whose ground state has energy gap as shown by Merrnin and Wagner. 28 
Thus, it is clear that the WZW term plays a crucial role to modify the ground state in 
a symmetric phase. Soon after, it was understood that the spin | chain is described by 
the 0(4) WZW-NLSM, and the onset of the valence-bond-solid order is understood by a 
marginal perturbation, which is interpreted as competition between spin and valence-bond- 
solid order. 29 Based on the 2D results, the WZW-NLSMs in higher dimensions are suggested 
to be in a different universality class from the LGW theory applying to various competing 
order physics.®ESH2Hll 

It is, however, significantly more difficult to analyze the models in higher dimensions 
than 2D partially due to a lack of a local conformal symmetry. Furthermore, a continuous 
symmetry can be spontaneously broken in higher dimensions, so a symmetry-broken state 
becomes another stable fixed point, which makes renormalization group flows complicated. 
The powerful theoretical tools such as e = D — 2 or large N methods are not applicable 
since changing space-time dimensions and an order parameter manifold are prohibited by 
the presence of the WZW terms. Therefore, it is quintessential to find ways to understand 
properties of symmetric phases in the WZW-NLSMs. Below, we find one concrete way to 
investigate symmetric phases of the WZW-NLSM type competing order physics at least in 
even-spacetime dimensions. 

Anomalies in WZW-NLSM 

We employ the standard strategy to investigate quantum anomalies: to promote a global 
symmetry to a local one (gauging) and search for inconsistency? 2 The promotion is done 
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by introducing a gauge potential and associate minimal coupling. Without the WZW term, 
the minimal coupling is enough to gauge the symmetry and there is no ambiguity. But, 
gauging the WZW term is subtle since it is expressed in X D+1 . One criterion of a proper 
gauging procedure is that gauging does not change space-time dimensions of dynamics. For 
example, equations of motions should be well-defined in the original dimensions (Xu) even 
after gauging since infinitesimally weak gauge coupling is conceivable. 45 

For notation convenience, we use differential forms following Nakahara’s^ (also, see Sup¬ 
plementary Materials). In 2D, the WZW term is 

Swzw = icz [ u 3 , u 3 = 4> % d(jA A d(j) k A dcp 1 . 

Jx 3 

Note that even though the WZW term is expressed in X 3 , the equation of motion is well- 
defined in 2D, 

1 i /r 

- + ^ K £tlU£ ijki^d ll (j) k d u (j) 1 = 0 . ( 2 ) 

This is because lu is the highest form (du = 0) in X 3 , and the closeness (du = 0) is tied to 
a two dimensional equation of motion from Poincare’s lemma. From now on, we focus on 
SO(D + 2) symmetry instead of 0(D + 2) since the’t Hooft anomaly matching only cares 
continuous groups. The variation of the order parameter in the SO (4) vector representation 
is Scj) 1 = e a [ta\)(tS with the Lie algebra [t„] t j = — [t a \ji of .SO^), (a = 1, • • • ,6). Local 
dependence of e a (x ) determines whether the transformation is global or local. 

We introduce a gauge potential A = A^tadx^ to gauge the SO(D + 2) symmetry whose 
held strength is F a = dA a + \ f^ c A b A A c . The covarinat derivative is 

Dcf) 1 = (d + A)(j) z = D^dx^ = (d^ij + dx*. 

Then, the minimally coupled WZW term in 2D is 

Swzw = ic 3 f u 3 , u 3 = 0* D(jA A D(j) k A Dcf) 1 . 

Jx 3 d! 

The Q notation is for the minimal coupling with a covariant derivative (12(5) —>■ O = D(D)). 
Gauge transformations (5 g ) of the potential and the covariant derivative are d g A a = — de a + 
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f bc A b e c and 5g{D(f) 1 ) = e a (x)[t a ]ij(D(fA). By construction, the minimally coupled WZW term 
is gauge-invariant. 

Its equation of motion from the minimal coupling, however, is not two dimensional. This 
is easily shown by an exterior derivative of the minimally coupled WZW form and we find 


du) 3 


-d 


F a A v, 


+ d ab (Q 3 ) ab 


^0. 


An one-form v a = ^£ijki[t„] t j(p k d(j) 1 and the three dimensional Chern-Simon (CS) term, 
( Q 3 J — A a A (dA b + Jcd -) are introduced. The anomaly coefficient is 

r _ £aia 2 f>i&2 


Thus, the simple minimal coupling is not enough to make the WZW term gauge-invariant 
and the equation of motion two-dimensional. 

To make the equation of motion well-defined in 2D, one can rearrange and find the closed 
form. The wedge product (A) is implicit hereafter. 


d 


^3 + F a V a + d a b (Q 3 ) 


\ab 


0. 


Thus, the gauged total action with a two dimensional equation of motion should be 

Stot = So + Swzw + (* £ 3 )^ J (F a Va + d ab . 

Note that the last two terms contain gauge potential and field strength, so they vanish 
without gauge fields. 

The gauged action, however, has significant inconsistency under gauge transformations 

SgStot = ~(ic 3 )k [ d ab de a A b , (3) 

Jx 2 

whose origin is the presence of the Chern-Simon (CS) term. Therefore, it is impossible to 
gauge the 50(4) symmetry in the WZW-NLSM, which indicates the symmetry is anomalous. 
Note that similar inconsistency appears at the boundaries of the 17(1) CS theory in quantum 
hall systems.® 
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In 4 D with the SO(6) model, the similar procedure is applied with little more tedious 
calculation. The five dimensional volume form is 

u = ft dft dft dft dft 1 dft \ 

5! 

and the exterior derivative of the minimally coupled WZW term is 

d(u + F a v a + F a F\ ab) + cU c Qf c ) = 0 


with the anomaly coefficient 


£ a\a,2b\b2C\C2 


The five dimensional CS term is (Qft abc = A a dA b dA c + ~ l f^A e A^A a dA b + 
AA a f b de A d A e fftA 9 A h ), and the three form is v a = F£ ijk i mn (t a )ift k dft dft deft. The sym¬ 
metrized interior derivative is introduced, i( a Vb) = di\ ab y It is straightforward to show the 
SO{ 6) symmetry is anomalous. 

Note that the anomaly coefficient is only non-zero when all the indices are used up. This 
indicates that gauging the full 5'0(4) in 2D and SO(6) in 4 D is crucial to find anomalies. 
Gauging subgroups such as SO (3) does not use up all the indices, so the coefficient automat¬ 
ically vanishes. In other words, gauging subgroups is always well-defined, and the anomaly 
structure only appears when the full symmetry is gauged. 

Three remarks follows. First, our anomaly coefficient calculation in WZW-NLSMs is 
independent of the coupling constant strength ( g 2 ). This is consistent with non-perturbative 
nature of quantum anomalies : in the extreme limit with a very large (bare) coupling 
constant, one can imagine a symmetric ground state. Still, the anomaly structure must 
be there, so there must be massless excitation to reproduce the anomaly structure since 
order parameters are completely energy-gapped in symmetric phases. The presence of 
massless excitation can be more rigorously shown by investigating singularity structure 
of currents correlation functions 47 ^ (see also Supplementary Materials). Second, the 
presence of massless excitation in 2 D, which is consistent with the Witten’s bosonization 
results, and 4 D with anomalies indicates their ground states are qualitatively different 
form LGW theory’s. Thus, our anomaly calculation shows the WZW-NLSMs describes 




massless symmetric phases qualitatively different from the LGW theory and provides 
non-perturbative criteria to distinguish quantum criticalities. Third, our calculation can 
also be applied to relations between topological boundary and bulk phases with non-abelian 
symmetries. The gauge-invariance and equation of motion properties are intrinsically 
connected through the presence of the D + 1 dimensional Chern-Simon terms. Their 
presence appears in a certain class of D+l dimensional topologically ordered phases. 
Therefore, our work explicitly shows that the D+l dimensional (non-abelian) topological 
phases, related with the non-abelian Chern-Simon theories, can have D dimensional gapless 
boundaries guaranteed by the presence of quantum anomalies. 

Anomaly matching and Minimal model 

Quantum anomalies guarantees the presence of massless excitation in symmetric phases 
but it does not pin down the symmetric ground state completely. But, at least, it is obvious 
that the conventional strong coupling analysis as in the LGW theory is not applicable to 
the competing order theories with the WZW term. 

A priori, it is not even clear how many phase transitions are and whether they are first or 
second orders in strong coupling limits. It is because all ignored higher order terms become 
important in strong coupling limits. Thus symmetric ground states of the WZW-NLSMs 
are not uniquely determined without further microscopic information. But, no matter what 
happens, the presence of massless excitation is guaranteed by quantum anomalies. 

We first consider one specific model which reproduces a given quantum anomaly structure 
which helps us to understand qualitative differences better. Then, we use’t Hooft matching 
condition to investigate candidates of symmetric phases. The minimal model to describe a 
symmetric phase with quantum anomaly can be obtained from the symmetry breaking pat¬ 
tern in the symmetry-broken phases, SO(6)/SO(5) ~ S 5 in 4 D. It is well-known in literature 
that the two color (N c = 2) QCD with the two flavor (Nf = 2) enjoys the enlarged -SC/(4) 
flavor symmetry instead of SU(2)l x SU(2)r. 51 53 The symmetry is spontaneously broken 
to Sp( 4) by the chiral condensate. Since S77(4) is isomorphic to SO( 6) and Sp( 4) to SO( 5) 
in terms of the Lie algebra, the symmetry breaking pattern is exactly SO(6)/SO(5) and the 
dynamics of the Goldstone boson is described by a SO (6) NLSM (see also Supplementary 
Materials). In 2D, it is well known that spin 1/2 chains realize the *S'0(4) WZW-NLSM. 

Inspired by the hints, we construct fermion models coupled to the SO(D + 2) bosons. 
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We introduce a complex spinor, T, which couples to the order parameters, 

Sf = J V\d T ~i'fds)V + (4) 

s is for spatial dimensions, and i is for order parameters. ( 7 s , Tj) matrices satisfy the Clifford 
algebra, s — 1 (1, 2, 3), and i — 1, • • ■ , 4 (1 • • • 6 ) in 2D (4D). By the Yukawa coupling, the 
fermions become massive in the symmetry broken phases. The minimum numbers of spinor 
components for the Clifford algebra are four in 2 D and sixteen in 4 D as in the spin chain 
and N c = 2 QCD. 

The WZW term can be easily reproduced in the symmetry broken phases by integrating 
out the fermions, 

T eff = — log Z, Z[4>] = f exp (—Sf) = DetD 

J 'll), 

with V = d T — i^ s d s + A0Tj. Note that matching the level of the WZW term automatically 
satisfies the anomaly matching condition, and we indeed find the same level by integrating 
out the fermions. In Supplementary Materials, standard field theoretic consideration with 
group structures 1 2 is presented to be self-contained which is useful in symmetric phases. It 
is also straightforward to show further gradient expansion in the symmetry broken phase of 
the minimal model reproduces the WZW-NLSM. 

The minimal model in 4 D is naturally constructed by fermions and bosons, 

s min = j <t\d T -i^d r )<t + 

+ (5) 

T is in 50(6) spinor representation (sixteen complex fields) and r is a tuning parameter. 
In 2 D, a similar model can be constructed which is well understood, so we focus on 4 D 
from now on. By tuning r, one can access the symmetry broken phase (r < r c ) and the 
symmetric phase (r > r c ) where r c is a critical value and its numerical value depends on a 
renormalization group scheme. The minimal model is a class of the so-called Higgs-Yukawa 
theory, and uniqueness of the minimal model is the symmetry structure of bosons and 
fermions associated with quantum anomalies. 
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FIG. 2: Schematic RG flow of the minimal model in 4D with two parameters. The horizontal 
(vertical) axis is to characterize fluctuation (anisotropy) strength. N + M is SO(6 ) symmetry 
broken phase characterized by 5 Goldstone modes. M (N) is SO(M ) ( SO(N )) symmetry broken 
phase characterized by M — 1 (IV — 1) Goldstone modes. S is SO(6) symmetric phase with massless 
excitation enforced by quantum anomalies. 

The renormalization group flow of the general Higgs-Yukawa theory is well understood ,® 4 
and our minimal model has the same structure. The model is at the upper critical dimension 
4 D, so the mean-field description is valid with logarithmic corrections. There are three 
fixed points : the symmetry broken phase N + M, the quantum critical point QC, and 
the symmetric phase S (fermionic) as illustrated in the horizontal line of Fig. 2. Note 
that the symmetry broken phases in the minimal model and the LGW theory are similar, 
but quantum critical points and symmetric phases are fundamentally different due to the 
presence of fermions. 

We consider perturbations which break the full SO ( 6 ) symmetry down to its subgroup 
H = SO{N) x SO(M), with N + M — 6 to connect the anomaly structure with competing 
orders (we treat S'O(l) ~ Zt). For example, SO(3) x SO(3) symmetry in the minimal model 
allows anisotropy operators 

OT - (</> i+3 ) 2 , Y (^ tr ^) 2 - (T f r* + 3 T) 2 . 

i=l,2,3 i=1,2,3 
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Again, the minimal model is at the upper critical dimension, so one can easily read off 
scaling dimensions of the operators at each fixed point. Near QC , the former operator is 
relevant because it reduces the number of massless modes but it is irrelevant at S since the 
order parameters are gapped. The four point fermion interaction is irrelevant in all three 
fixed points. Schematic RG flow of the minimal model with the anisotropy parameter is 
straightforwardly obtaind as in Fig. 2. Different symmetry breaking patterns with different 
N, M give similar RG flows. Note that the RG flow structure of the minimal model is 
similar to one of the LGW theory, but crucial distinction between the minimal model and 
the LGW theory is the presence of massless fermions in S enforced by quantum anomalies. 
If the fermions are identified as electrons, the symmetric phase is nothing but Weyl or Dirac 
semi-metals. 5 

Before closing this section, we emphasize that the reason we consider the minimal model 
is to provide a concrete example of the symmetric phases. But, our anomaly calculation is 
powerful enough to be applied to more generic cases even with strongly correlated ground 
states. We discuss such more generic cases below. 

Non-minimal models 

Let us consider non-minimal models which have the same anomaly structure. A pri¬ 
ori, all conformal held theories (CFTs) with the same anomaly coefficient are candidate 
theories of S. One straightforward way to construct non-minimal models is to consider dif¬ 
ferent representations of 577(4) fermions instead of a single 5 , D(4) fundamental representa¬ 
tion of the minimal model. Detailed discussion about other representations is presented in 
Supplementary Materials. Notice that the RG how structure of the models with different 
representations is basically the same as the minimal model’s. 

If fermions are not electrons but fractionalized particles such as spinons, then the sym¬ 
metric hxed point can be identified as a spin-liquid phase. If spinons are weakly coupled to 
gauge held such as 1/(1) or Z 2 , then all the properties of the electronic minimal model is 
inherited and the symmetry broken phases contain remaining gauge structure, so-called * 
phases. They do not describe conventional symmetry broken phases. Therefore, if the sym¬ 
metric hxed point with spinons are adjacent to conventional (confined) symmetry broken 
phases, the gauge structure must be non-abelian. For example, the well-known Banks- 
Zak hxed point— with spinons could be a candidate of the symmetric hxed point S. Then, 
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FIG. 3: Marginal and relevant operators with H = SO(N) x SO(M) symmetry, (a) RG flow 
with a marginal operator around the symmetric phase (top). The horizontal (vertical) axis is to 
characterize fluctuation (anisotropy) strength. Phase diagram with a marginal operator (bottom). 
Sign of the tuning parameter determines relevance (or irrelevance). Thus, one side of phase diagram 
is a 50(6) symmetric CFT, and the other side (A") is either a H symmetric CFT or a H symmetry 
broken phase, (b) RG flow with a relevant operator around the symmetric phase. The relevant 
operator drives the RG flow from 5 to Xm and X at. They are either H symmetric CFTs or 
H symmetry broken phases. Once Xm and X jy break different symmetries, the symmetric fixed 
point (5) may describe a “deconfined” quantum criticality. The curved (black) lines in the bottom 
figures are for energy scales associated with symmetry breaking (e.g. order parameter scale). 

condensing the order parameter endows spinon energy gap, and the remaining gauge field 
becomes confined naturally. 

In principle, two different RG flow structures near S are possible if the symmetric fixed 
point S is described by strongly coupled CFTs with the same anomalies. In contrast to the 
minimal model’s RG flow, anisotropy operators could be relevant or marginally relevant. 
If marginally relevant, then the symmetric phase becomes stable with one definite sign of 
the coupling constant, but the opposite sign makes the symmetric phase unstable. Thus, 
one side of phase diagram is described by a 50(6) symmetric CFT, and the other side (A") 
is described by either a H symmetric CFT or a H symmetry broken phase as shown in 
FIG. [3](a). If relevant, then the symmetric phase becomes unstable with both signs of the 
coupling constant. Again, the final states can be either //-symmetric CFTs or //-broken 
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phases as shown in FIG. [3](b). Once the final two states break different symmetries, then the 
symmetric fixed point (S) connects two broken phases directly, which describes deconfined 
quantum criticality® 

The above discussion gives the two necessary conditions to realize deconfined quantum 
criticality in SO(6 ) WZW-NLSM in AD : anomalies and relevant symmetry breaking 
operators. These conditions provide non-perturbative criteria to characterize universality 
class of quantum phase transtions between competing orders. For example, the universality 
class of the SO (Jo) WZW-NLSM cannot be the same as one of the non-compact CP (A) 
model in AD due to the absence of anomalies in the latter model. Notice that in 3 D the 
CP (A) model and the 50(5) WZW-NLSM are proposed to describe the same universality 
class but our anomaly criteria do not applied to odd space-time dimensions. 

Discussion and Conclusion 

In experiments, direct measurement of the non-abelian anomalies associated with com¬ 
peting orders is even more difficult than one of the chiral U (1) anomaly in Weyl semi-metals 
because we do not know how to couple the non-abelian current directly in experiments.® ® 
Yet, there are traits associated with the anomalies. 

Protection of massless excitation is one of the most significant characteristics of the 
presence of quantum anomalies. Their numbers are, however, not universal. For example, 
in the minimal model in AD, the numbers of massless excitation in Goldstone phase, quantum 
critical point, and symmetric phase are 5 (bosons), 21 (= 6 + 16) (bosons +fermions), and 16 
(fermions). In non-minimal models, the symmetry broken phase has the same number, but 
critical point and symmetric phase have different numbers of massless excitation. Clearly, 
this is different from LGW theory’s where all massless excitation has definite numbers at 
each fixed point. The different numbers of massless excitations may contribute to transport 
differently, which is in principle measurable. 

In the minimal model, characteristics of anomaly becomes more evident. First of all, 
semi-metallic behaviors appear in a symmetry restored phase or high temperature regime 
(but lower than cut-off scale, say band width) if fermions are electrons. Massless elec¬ 
trons in a symmetric phase (or quantum critical regime) governs low energy physics, so 
electrical and thermal currents are carried by electrons with the linear spectrum. Natu¬ 
rally, the Wiedemann-Franz law holds especially in non-hydrodyanmic limits. By lowering 
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temperature, the £77(4) symmetry can be broken and the electrons become gapped. Thus, 
electrically, insulating behaviors (energy gap) appear, and the symmetry breaking transition 
is concomitant with the transition between semi-metal and insulator. On the other hand, 
Goldstone modes from spontaneous symmetry breaking carry thermal currents even though 
electrons are gapped. Since Goldstone modes and massless electrons have same dispersion 
relation, thermal transport in symmetry broken phases is qualitatively similar to the one in 
symmetric phases. Thus, near the semi-metal and insulator transition, electric and thermal 
currents behave differently, and the Wiedemann-Franz law would be violated. 

One experimentally realizable system of the Si7(4) anomaly is pyrochlore systems with 
all-in-all-out magnetic order parameter 59 in addition to N c = Nf = 2 QCD . A class of 
pyrochlore structure is described by a quadratic band touching modef 60 and the onset of 
all-in-all-out order parameter 61 induces eight Weyl points (16 fermions), the minimal nec¬ 
essary number to realize the SO(6) ~ SU(4) anomaly. We note that evidence for all-in 
all-out ordering and violation of the Wiedemann-Franz law in spin-orbit coupled pyrochlore 
structures has been reported in literature 62 63 though precise connection with anomalies need 
further investigation. 

In this paper, we investigate non-abelian anomalies in quantum phase transitions with 
competing orders. We show that the WZW-NLSMs in 2D and 4 D have quantum anomalies 
by calculating the anomaly coefficients. Non-perturbative nature of the anomalies allows 
us to investigate not only a symmetry broken phase in weak coupling limit but also a sym¬ 
metric phase in strong coupling limit even though the presence of the WZW term prohibits 
conventional e = D — 2 and large N expansion methods. Applying the ’t Hooft match¬ 
ing condition, it is shown that the universality class of the models is qualitatively different 
from the conventional (p A theory’s. In sharp contrast to the gA theory, symmetric ground 
states of WZW-NLSMs contain massless excitation though their numbers are not uniquely 
determined. Thus, we construct the minimal model and investigate its properties under 
anisotropy perturbation. Then, we extend the model to more general ones and discuss im¬ 
plication of the anomalies in competing order physics. Further research on anomalies, for 
example, parity anomaly in odd space-time dimensions and mixed anomalies with gravity 
(thermal effects) in connection with topological phases (NLSM with the theta term) are 
desirable. 
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Appendix A: Massless excitation with anomalies 

It is well understood that massless excitation is guaranteed by continuous symmetry 
anomalies)- 4 ' 49 The presence of continuous group’s anomalies enforces singularities of ana¬ 
lytical structures of currents correlation functions. To be self-contained, we introduce the 
proof with slight modification following the notation in Coleman and Grossman . 49 
In AD, the anomalous Ward identity is in three currents correlation function, 

iVA0i,g2,g3)£ (4) (<?i + <72 + <73) = 

[ < 0\Jfi{Xl)J v (x2)Jx(xg,)\0 >, 

i 

and the current conservation gives 

Qs) = Ae^ af} q*q%. (Al) 


All non-abclian Lie algebra indices are absorbed into the anomaly coefficient A. 

The correlation function is symmetric under simlutaneous permutations of (qi,q 2 , < 73 ) 
and (/ 1 , v, A). Now let us investigate analytic structure of the correlation function. Due to 
permutation and covariance, the structure must be in the form 


IVa = F(q 2 i ) 


e#«/a09l 92 03A + ^Aa/^IiW + W/30301 Q2v 


We omit possible tensors which cannot contribute to the anomalies. Note that the mornen- 
tums are off-shell, so one can access all available regions and we focus on the region 


2 2 2 
01 = 02 = 03 


-Q 2 . 
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The correlation function contracted with q%\ gives 


?3 r ^A(gi, 92, 93 ) = ~F{Q 2 )Q 2 e flua ^q^. (A2) 

Then, the anomaly equation (A.l) gives 

OQ 2 ) = 

The pole structure at zero mass nicely show the presence of massless excitation (also sed^ 8 
for dispersion analysis). The singularity even further enforces that UV and IR information 
needs to be matched. 

In the paper by Coleman and Grossman, they add more conditions such as non- 
singularties from vertex corrections, and they conclude the helicity of massless degrees of 
freedom is ±|, which indicates the symmetric phase is massless fermions as in our minimal 
model. The authors argue that the assumptions are not that strong, so it would be very 
interesting the conditions are proved / disproved in future research. 

The above discussion only relies on the anomaly properties and nothing more, thus it is 
applied to everywhere in phase diagrams. But, it is only applied to anomalies of continuous 
symmetries since the current conservation plays a crucial role. For the discrete gauge group, 
which is especially important in SPT physics, the presence of anomalies does not guarantee 
massless excitation. 8 11 

We note that in 2 D, the minimal symmetry for spin 1/2 chains to be massless is SU (2) x Z -2 
corresponding SU(2 ) x Z 2 16 which is smaller than 50(4) ~ 577(2) x 51/(2), and it is manifest 
some subgroups of the continuous group is enough on lattice systems, and it would be 
interesting to find criteria to determine the subgroups in higher dimensions. 

Appendix B: N c = 2 QCD theory 

In QCD, confinement issue is subtle and it is known that chiral femrions with non-abelian 
gauge fields are confined in the infra-red limit. To connect the QCD theory to our minimal 
model, one needs to control one parameter to access both the chiral symmetry broken phase 
and the symmetric phase. One way to do this is to tune gluon-quarks interaction strength 
by adding a color doublet, flavor singlet scalar (Higgs) field < F. In one limit, $ is massive 
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and then the chiral symmetry is broken by confinement, but in the other limit, the Higgs 
field obtain vacuum expectation value (vev) and the gluon fields become screened. A phase 
transition between the two limits can be reached by dialing the Higgs vev. Due to SU( 4) 
anomalies, one now expects physical massless particles in the spectrum. They are described 
by colorless product of the quark field and the Higgs. Thus, in the symmetric phase, one 
has massless fermions in the fundamental representation of 577(4), which could be identified 
by our minimal model. 

Below, we provide more information on symmetry enlargement and order parameter 
construction in N c — Nf — 2 QCD suppressing gluon fluctuations. 

The massless Dirac Hamiltonian of N c = Nf = 2 QCD in the chiral representation is, 

'Ho = q\i o k d k q L - q ] R i v k d k q R . (k = 1,2, 3) (Bl) 

q rR is a spinor with Nf x N c x 2 components (flavor x color x Lorentz). The flavor indices 
are implicit and it is manifest that the Hamiltonian has SU^Nf) x SU R (Nf). 

For N c = 2 QCD, the symmetry is enhanced to a bigger symmetry, SU(2Nf). One can 
see this by employing the chiral particle-hole transformation for a chiral field (say, q R ), 

q R = {T 2 a 2 ){q ] ) T ', q ] R = -q T (T 2 a 2 ), q = q L . (B2) 

Ti 2,3 are anti-hermitian generators for the 577(2) color. We define the corresponding Her- 
rnitian operators, t k = iT k . In terms of a new fields T 7 = (q, q), the Hamiltonian is 

H 0 = q\ i a k d k q L - q ] R i cr k d k q R 

= q ] i cr k d k q + qU a k d k q = T 1- ia k d k T. 

SU(2Nf ) symmetry is manifest with new spinors. Notice that all fermions are chiral in this 
representation. 
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One mass term (color singlet and flavor singlet) in two different representations is 


QrQl + vUr = {T 2 a 2 )q + q\T 2 o 2 )(q ] ) T 

= -i q T (t 2 a 2 )q - q ] {t 2 a 2 ){q ] ) T 
= '& T p 2 t 2 a 2 '& + ^p 2 t 2 cr 2 (^) T 


(B3) 


More mass terms can be systematically obtained by using the Majorana spinors, 


'I' - 'll - i 'll- '! = ('/1 'h! 


(B4) 


Then, the fermion bilinears have the form 




T/ M-M t n M + M t 
V (-g^- 2 - 

p T M{p 3 + ip 1 )'?] 


= 77 T M(/i 3 — ip 1 )^ 

7fMp 3 p = + T t M(T t ) T ) 

q T Mp l, q = — (T r MT - 


with /z 0,1 ’ 2 ’ 3 are for the Majorana space. Thus, the mass term in the Majorana representation 

is 


QrQl + q[qR = ^ p 2 t 2 a 2 (^) T ) 

= rj T (t 2 a 2 pV)r). 


The kinetic Hamiltonian with the Majorana representation is 

Uq = p T i'y k d k rj 

= q T i(a 1 d 1 + a 2 p 2 d 2 + a 3 d 3 )rj. (B5) 
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For the reference, we summarize the Pauli matrix symbols, 


f : color, a 1 : Lorentz, p 1 : (q, q), 

p l : Majorana, r l : flavor (B6) 

and the matrices of the kinetic term are 

7 1 = t 0 a 1 p°n°T 0 7 2 = t°ij 2 p° p 2 r° 7 1 = tVp° A° 

All zero components of the matrices are the identity matrices in the corresponding spinor 
spaces. 

Three more mass terms can be obtained in the flavor triplet channel. 

* - * Qr^Ql, (j = 1, 2,3). (B7) 

In the Majorana representation, the mass terms are 

T/j.2 2 2 1 1\ T( ±2 2 1 3 2\ T14.2 2 2 1 3\ 

rj (t a pprjrji rj (t a p p t ) q, rj (t a p p r )rj. 

Two more mass terms are from superconducting channels 

{q T L t 2( r 2 T 2 q L + q[t 2 a 2 r 2 (q[) T ) + (l -7 R^j, 
i(q T L t 2 a 2 T 2 q L - q[t 2 a 2 t 2 (q[) T ^j - (l ->• r) , 

and in the Majorana representation they are 

q T (t 2 a 2 p 3 p 3 r 2 ), rj T (t 2 a 2 p° 

Thus, we have six mass terms, rj T V a r] with a = 1, • • • ,6, 

tU _ + 2 2 2, 3 0 ti 2 +2 2 2 1 1 r 3 _ +2 2 1, 3 2 

1 — t (J p fJ, T ) 1 — t (7 p p T , 1 — t (J p pj T j 

t^4 4.2 2 2 1 3 -n5 4.2 2 3 3 2 -p6 ,2 2 0 1 2 

I = t a p p t , 1 = t a ppr, T — t a p p r . 
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Note that all gamma matrices ( 7 k , T a ) anti-commute each other (Clifford algebra) and all 
six order parameters are color-singlet. Thus, our construction is independent of introducing 
the Higg’s doublet to control the gauge coupling strength. 

The mass term Hamiltonian is parameterized by six fields, 

'Hm — gr] T (f) a r a r], (B 8 ) 

with a dimensionful coupling constant g. The six fields are matched with the 517(4) gener¬ 
ators from which the symmetry breaking down to SP( 4) is described by 

HI - s 5 -► rr = 1 . (B9) 

The last condition guarantees the correct ground state manifold. 

Appendix C: Anomaly coefficients and non-minimal models 

One useful diagnosis of the presence of anomalies is provided by ’t Hooft . 3 Anomaly 
appears when a global symmetry is gauged and its gauge transformation produces physical 
effects, which are inconsistent with definition of gauge transformations. Such inconsistency 
with gauge transformations must be canceled because infinitesimally weak gauge coupling 
is always conceivable. To cancel it, one can add non-interacting massless fermions (so- 
called spectators) with opposite anomalies. Then, after spontaneous symmetry breaking or 
confinement, the original sector should contain some massless degrees of freedom to cancel 
the anomalies of the non-interacting fermions. 

Below, we follow the previous discussion in 577(4) representation to consider the anoma¬ 
lies in different irreducible representations . 7 A representation R of the SU(4) group is char¬ 
acterized by three numbers qi, q- 2 , < 73 , which are the numbers of columns in the Young tableau 
with 1, 2, and 3 rows. 

The dimension of the representation is 

D{R)= -^(?i + 1)(? 2 + 1)(<Z3 + 1) 

x {Qi + 92 + 2) (<72 + <73 + 2) (qi + <72 + q3 + 3) 
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The anomaly is characterized by the number A(R), 


Tr({T“, T b R }T c R ) = A(R) Tr({f“, t b }t c ) (Cl) 

where t a refers to the fundamental representation of SU(4). 

A(R) = ^r(gi - qs)(qi + gs + 2)(gi + 2 q 2 + q 3 + 4 )D(R) 

60 

If R — (qi, q 2 , qs) then R = ( q 3 , q 2 , qi) and A(R) = — A(R ). For real representations q\ = qs 
and A(R) = 0. In the table below we exclude the complex conjugate representations. 


with d([ = 4 and A ^1 = 1. 


) can be constructed by two conditions 
: 1) anomaly matching and 2) presence of Yukawa coupling (fermion mass in the symmetry 
broken phase). The formal expression for the anomaly condition is 

J^n R A(R) = k, (C2) 

R 

k is the WZW level and the summation is over all representations. Possibility of the Yukawa 
coupling can be obtained by multiplication of representations. 

Appendix D: Differential Geometry 

The interior derivative is defined as 

i.fi = (Dl) 



A class of non minimal models (weakly coupled one 
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with a vector field £ a = (j) z {t a )ijdj. 

Since the volume form is the highest form, we have nice properties 


du = 0, i a uj = dv a , C a {uj ) = 0 


We introduce 0 notation for replacing a partial derivative with a covariant derivative 

(Q(d) -> £l(D)). 

For 1 + 1 dimensions, one can show that 

*o(w) = ^e i jki(j) l i a (d(j> ] )d(j) k d(j) 1 -\ - 

2 21^ijki{ta)ijd(fr d<j> dv a 

which defined the one form, 

+a g d(j) . 

The relation (t a )ij = \eijki{u a )ki is useful. One can show C a Vb = fab v c straightforwardly. 
The interior derivative of the one form is 


h ( va ) = 2 ((/‘) 2 + (/ 2 ) 2 )- 


(D2) 


Its symmetrized one is 


iqjvb) + ib(v a ) 
2 


1 

4 


^aia2?>i&2 


One useful identity is 


(D3) 


,l<! = ,!<)- 


One can show this by using components of the form. By using the identity twices, one can 
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obtain 


du = —F a dv a = -F a (dv a + F b i b v a - A b C b v a ) 
= -d{F a v a ) - F a F\ b v a) . 


Thus, the anomaly coefficient is 


, _ i a (vb) + ib(v a ) _ 1 

® ab 2 ^<110261 &2 


(D4) 


For 3 + 1 dimensions, one can show that 


v a 2 ^| '- ijkhnni^a)ij P dtp dp dp , C, a V b fab^c 


By using the identity, we obtain 


du = -d(F a v a ) - F a F b d^T b) 

= -d{F a v a + F a F b V(Pp) - d abc F a F b F c . 

The relation i^ a v b ) = di\ ab ) is used, and the anomaly coefficient is 


dabc 


*c^(o6) + ibV( oc) i'a'V(cb) £010261620^2 

3 = 4! 


(D5) 


Appendix E: WZW term from integrating out fermions 


Main purpose of this section is to provide explicit informaiton on fermion determinant 
calculation. 

In 1 + 1 dimensions, let us consider the four component complex spinors, T, with two 
pauli matrices The action is 


S 


I ¥{d T -iaa z T°d x )ty + 


with anti-commuting matrices {M,, Mj} = 25ij and {a z ,Mj} = 0. The coupling constant 
a = ±1 will determines the sign of the WZW term. In the next senction, the mass matrices 



(Mi) are explicitly constructed. Since the matrix size is 4 x 4, the total number of the 
matrices is five. 

The effective aciton is defined 

T eff = - log Z, Z[<t>\ = I exp (~S f ) = DetD 

with V = d T — iaa z T°d x + gcj) and (j) = <p l M t . 

The variation of the effective action is 

-5^Y e 2 ff = tr (5DD\DD^)~ l ) 

= g 2 tr(<50 ■ cj)Go(l + MGq + •••)). 

The Green’s funciton is introduced by the operator, 

DD t = — d 2 — d 2 + g 2 + g(d T (fi — iacr z d x (f)) = Gq 1 — M. 

Gq l (k,ui) = u 2 + k 2 + g 2 . With this set-up, we do the gradient expansion in terms of g (or 
M). 

In 1 + 1 dimensions, the topological term can be obtained by the second term 

-V V { ? = g 2 tr( 5 j) ■ j)M 2 G 3 0 ) 

= i [ ^-£^e ij ki5(j) l cty i d ll (t) k d u (j) 1 . 

Jx 271 

The relation tY(M l MjM k Mia z ) = 4 e^i is used. From this calculation, one can further 
notice that the sign of the imaginery term depends on the sign structure of the Hamiltonian. 
Instead of Hq = 4d(— ia z r°d x )^ 

In 3 +1 dimensions, let us consider the 16 component complex spinors, T, with four pauli 
matrices. There are nine gamma matrices for the Clifford algebra which are construced in 
the previous section. The N c = 2 QCD with Nf = 2 theory has the same number of degrees 
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of freedom. Following the previous calculation, we introduce 


S = j ¥ (d T - ia Yd s )^ + gtf&riV 
-T eff = log D — log(<9 T — iarfd r + gfi), = 0T ?; 

with s — 1 , 2 , 3 for the gamma matrices in the kinetic term and i, j = 1 , ■ ■ ■ 6 for the gamma 
matrices in the mass terms ({r^Tj} = 26^). The variation of the effective action is 

_ S^T eff = tr (SDD'(DD')- 1 ) 

= g 2 tr(8(f) ■ 0Go(l + MG 0 + •••)) 


with 


DD 1 = — d 2 — d 2 + g 2 + g(d T (j> — iaj r d r $) = G 0 1 — M. 

In 3 + 1 dimensions, the topological term can be obtained by the fourth term 
- 5^ r< 4 > = g 2 tr(5j) ■ 0 M 4 G 5 0 ) 

, which is the exactly same as the variation of the WZW term. 

Appendix F: Group theory consideration of anomalies 

In this section, we provide more information on group theoretical relations in non-abelian 
anomalies. The group theoretical relation is especially powerful to find anomalies in fermion 
only theories. One can calculate either a loop-diagram or the Fujikawa measure in path 
integral. 1 2 

1. 1+1 dim : SO( 4) ~ SU(2) L x SU(2) R 

We use the Clifford algebra to connect representations of S'0(4) group and SU(2) L x 
SU(2)r. The idea is SO( 5) group has a natural spinor representation in terms of five 
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Gamma matrices, and we use four gamma matrices out of the five. The remaining one 
becomes the “chiral” operator. To construct it, we use 

r 5 = AV 0 

T^aV, r 2 = a 2 r 2 , T 3 = cr 2 r 3 , T 4 = aV 0 

It is straightforward to relate them to the Lie-algebra of 50(4). We define M tJ = |TT J . 
Explicitly, we have 

M 12 = a \ 3 , M 13 = l -a\ 2 , 

M u = l -a*T l = T 5 M 23 , M 24 

M 34 = l -a 3 r l = T 5 M 12 

Note that T 5 indeed plays a role as the chiral operator. 

The anomaly coefficient is related to 

d ab = Tr(T 5 M a M 6 ) = £ aia2blb2 . 

2. 3+1 dim : 50(6) ~ 50(4) 

As in 1 + 1 dimensions, we can use the 50(7) representations. Seven gamma matrices 
are necessary, and their minimum size is 8 x 8. One needs three types of Pauli matrices and 
we use the notation (<t q/37 = a a ® cO <E> cr 7 ) with Pauli matrices cr a , a = 0,1, 2, 3. cr° is the 
2x2 identity matrix. Similarly, (a a/3 = a a <E> a 13 ). 

From the 50(7) algebras, we pick one operator as a projection operator. 

p7 = ^300 

p6 = r 5 = ^220, r 4 = ^230^ 

p 3 = (T 10 3 , r 2 =fT 102 ; r 4 = (7 101 > 


M 23 = Vr 1 

= l -a \ 2 = T 5 M 13 
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With Mij = PP, we have 

l\/[ _ *003 /\/t _ * 002 j\/[ _ *331 /\/f _ *321 

1V1 12 — VJ , M 13 — IG , iK/i4 — 2(T , iK/15 — %G , 

M 16 = «t 311 , M 23 = ia 001 , M 24 = kt 332 , M 25 = ia 322 , 

M 26 = fa 312 , M 34 = ^ 333 , M 35 = ^a 323 , M 36 = m 313 , 

j\/[ _ *010 i\/f _ *020 /\/r _ *030 

1V1 45 — Z(J , -^46 — 2(T 5 ^56 — ^(T 

Note that T 7 indeed plays the chiral operator. Thus, we can choose either +1 or —1 com¬ 
ponent and obtain the spinor representation of 50(6). 

M \ 2 = a 03 , M 13 = a 02 , Mu = a 31 , M 15 = a 21 , 

Miq = a 11 , M 23 = a 01 , M 24 = a 32 , M 25 = a 22 

M 26 = a 12 , M 34 = a 33 , M 35 = a 23 , M 36 = a 13 

M 45 = a 10 , M 46 = a 20 , M 56 = a 30 . 

These are also SU (4) fundamental representations. The anomaly coefficient is proportional 
to 

dabc Tr(.\/„{.!/,, A /,.} ) £ a\a,2b\b2C\C'2 

upto a non-zero numerical constant. 
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